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PART  1 


ELECTRON-IMPACT  IONIZATION  OF  THE  OXYGEN  ATOM 
1.  INTRODUCTION 

Electron-impact  ionization  of  the  oxygen  atom  is  an  important  basic  process  in  at¬ 
mospheric  and  plasma  physics.  Experimental  measurements  of  ionization  of  atoms  and 
molecules  including  oxygen  atom  have  been  reviewed  by  Kieffer  and  Dunn  [1].  More  recent 
complications  of  ionization  data  are  given  by  Laher  and  Gilmore  [2],  and  by  Bell  et  al.  '.3  . 
The  electron-impact  ionization  cross  sections  of  the  oxygen  atom  were  measured  by  Fite  and 
Brackmann  [4],  Rothe  et  al.  [5],  Brook  et  al.,  [6]  and  Zipf  [7].  In  the  first  two  experiments 
[4,5]  a  partiaUy  dissociated  oxygen  beam  containing  O2  and  0  was  used  to  measure  the  ratio 
of  the  ionization  cross  sections  of  the  two  spedes,  <r(0'^)/o-(02'''),  which,  upon  combining 
with  the  total  molecular  ionization  cross  sections  of  Tate  and  Smith  [8],  gives  (t(O^).  Brook 
et  al.  [6]  used  the  charge-exchange  process  as  the  source  of  atomic  oxygen  beam.  Zipf  [7] 
used  a  beam  of  O/O2  mixture,  but  the  normalization  procedure  was  different.  Zipf  [7]  also 
revised  the  data  of  Fite  and  Brackmann  [4]  using  the  newer  O2  ionization  cross  sections  of 
Mark  [9]. 

An  early  calculation  of  ionization  cross  sections  of  oxygen  atoms  [10]  was  made  by  utilizing 
the  relationship  between  the  electron-impact  and  photoionizaton  cross  sections.  Subsequent 
calculations  [11-16]  involve  the  Bom  approximation,  and  in  all  but  one  case  [16]  the  hydro- 
genic  Coulomb  functions  were  used  to  describe  the  ionized  electrons  and  the  2p  atom' ;  wave 
functions  were  obtained  by  various  means.  The  Coulomb  functions  were  orthogcnalized  to 
the  2p  function  by  the  customary  (but  somewhat  arbitrary)  procedure.  Peach  [11]  used  the 
2p  function  of  dementi  [17],  and  computed  the  cross  sections  with  both  the  Bom  and  Born- 
Ochkur  approximations.  Other  kinds  of  2p  atomic  wave  functions  were  used  in  Refs.  13-15  as 
explained  in  the  respective  papers.  Burnett  and  Rountree  [16]  computed  the  ionization  cross 
sections  by  the  Born  approximation  with  no  allowance  for  the  electron  exchange,  but  with 
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a  set  of  wave  functions  much  more  elaborate  than  the  ones  used  in  other  works  [11-15].  For 
the  initial  state  of  the  target  atom,  a  six-term  configuration  interaction  (Cl)  wave  function  is 
used  to  describe  the  groimd  state  of  the  oxygen  atom  0(®P).  The  continuum  wave  functions 
for  the  ejected  electron  in  the  final  target  state  are  determined  by  solving  the  close-coupling 
equations  that  result  from  the  inclusion  of  the  '*S,  ^D,  and  terms  of  the  2p^  configuration 
of  the  free  O'*'  ion  in  the  Schroedinger  equation  for  the  target. 

in  this  paper  we  conduct  a  systematic  theoretical  study  for  the  electron-impact  ionization 
of  the  oxygen  atom  from  the  excited  states  as  well  as  from  the  ground  state.  We  start 
with  the  Bom  approximation  to  calculate  the  ionization  cross  section  of  the  ground-state 
0(^P)  oxygen  atoms,  and  treat  the  electron  exchange  using  the  procedure  suggested  by 
Ochkur  [18],  i.e.,  the  Bom-Ochkur  approximation.  To  improve  the  cross  sections  at  low 
energies,  we  replace  the  Bom-type  approximation  by  the  method  of  (exchange)  distorted 
waves  with  Ochkur’s  exchange.  The  calculation  is  extended  to  some  low-lying  excited  oxygen 
atoms  in  the  2p^(^S)no^o  states.  To  calculate  the  wave  fiinctions,  we  use  the  frozen-core 
approximation  by  which  we  assume  that  Is,  2s,  2p  functions  remain  imchanged  in  the  0(^P), 
0[(^S)no^o  ^Lo],  and  0'*'('‘S,*D,*P)  states.  The  wave  functions  of  the  excited  orbitals  rioio-, 
of  the  ejected  electron,  and  of  the  distorted  projectile  electron  are  computed  by  the  Hartree- 
Fock  procedure  while  those  of  Is,  2s,  2p,  are  taken  from  the  paper  of  dementi  and  Roetti 
[19]. 

In  Section  2  we  set  the  theoretical  framework  under  which  the  cross  sections  are  ccdcu- 
lated.  In  Section  3  we  present  the  results  and  compare  them  with  the  previous  theoretical 
calculations  and  with  the  experimental  measurements.  Section  4  concludes  the  paper. 

2.  THEORY 

A.  General  Description  and  the  Bom-Ochkur  Approximation 

Since  the  theory  of  electron-impact  ionization  of  atoms  has  been  treated  extensively 
by  Peterkop  [20],  Rudge  and  Seaton  [21],  and  Rudge  [22],  we  will  cite  the  relevant  results 
developed  in  those  papers  without  detailed  explanations.  Let  us  first  consider  the  simplest 
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case  of  e-H  ionizing  collision  where  an  incident  electron  of  momentum  hko  impinges  on  the 
hydrogen  atom  of  wave  function  <^0(^2)  resulting  in  a  scattered  electron  of  momentum  hki 
and  an  ionized  electron  of  momentum  hk2.  The  kinetic  energies  of  the  electrons  are  related 
as 

kl-2I  =  l^  +  ]4,  (1) 

where  I  is  the  ionization  potentizd  energy  in  au.  The  bound-state  function  <^o(0  expresed 
as 

=  T-'^Pr^MYurr.A^),  (2) 

and  the  wave  function  of  the  ionized  electron  <^{k2,r)  as 

(3) 

Here  we  used  symbol  £,  which  is  k2/2,  to  indicate  the  energy  of  the  ionized  electron.  As 
indicated  in  Eq.  (2)  we  use  Pnj(r)  to  designate  the  “reduced”  radial  wave  function  rRni(r). 
The  cross  sections  of  this  process  are  dictated  by  the  direct  and  exchange  amplitudes,  f(fci  ,^2) 
and  g(fci,fc2).  Discussions  of  f  and  g  can  be  foimd  in,  for  example,  the  paper  by  Rudge  and 
Seaton  [21]  where  the  symbols  k  and  k*  are  used  in  place  of  our  fci  and  k2-  The  ionization 
cross  section  integrated  over  all  directions  of  ki  and  k2  is  given  by 

j  ={hk2lK)  j \f{kx,k2)±g  {ki,k2)\^dkidk2,  (4) 

where  the  superscripts  “-h”  and  on  the  left-hand  side  correspond  to  scattering  states 
with  total  spin  0  and  1  respectively.  If  the  H  atom  is  ionized  by  unpolarized  electrons,  then 
the  cross  section  is  the  weighted  average  of  the  “-I-”  and  with  the  statistical  weight  1:3. 

In  the  Bom  approximation  (with  no  exchange),  the  g  term  is  neglected  and  the  f  term  is 
approximated  by 

/(fcifcj)  ~  -(l/27r) y  <f>o{f2)  In  -r2|"*  4>"{k2,r2)  dfidfz,  (5) 


/  da 
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where  and  axe  defined  in  Eq8.(2)  and  (3).  Alternatively,  one  can  adopt  the 

Born-Oppenheimer  approximation  in  which  g  is  taken  from  Eq.(5)  by  permuting  the  two 
electrons  in  the  final  state,  i.e., 

g{kiX)  =  -(l/S^r)  J  4>o{^2)\ri  -  <^*(^2.^1)  (6) 

The  use  of  the  Born-Oppenheimer  exchange  amplitude  of  Eq.(6)  leads  to  the  well  known 
difficulty  of  overestimating  the  cross  sections  especially  at  incident  energies  near  the  thresh¬ 
old.  Ochkur  [18]  has  suggested  an  approximation  for  treating  the  exchange  amplitude  for 
ionization  processes,  i.e.. 


This  Bom-Ochkur  (BO)  approximation  is  shown  to  yield  much  closer  agreement  with  ex¬ 
periment  for  e-H  ionization  than  the  Bom  approximation  or  the  Bom-Oppenheimer  ap¬ 
proximation  [18].  In  the  present  work  we  shall  use  Eq.(7)  to  compute  exchange  amplitude 
yielding 

=l,k,k,/K)J(X*fmk,i^)\’dkdk.  (8) 

with 

=  l±\ko-k\V{kl-f4)-  (9) 

Following  Peterkop  [20]  we  evaluate  Eq.(8)  by  re-expressing  f(fci,fc2)  of  Eq.(5)  as 


fihM)  =  -2^--*  <  Oj^lfco  >, 

K  =  ito  —  kii 

<  0\K\ik  >  =  /  <^o{^2)e‘^  <f>*{k2yr2)df2. 


Integration  over 


dfci  and  dkj  reduces  Eq.  (8)  to 


<  njo  \K\€  >  \^K-^dK, 


(10) 

(H) 

(12) 


(13) 
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I  <  nolo\K\e  >\^  =  {2L  +  J  Ar)e'^'^^R.i{r)Yr^if)di^\  (14) 

The  integration  limits  Kmi„  and  K _ in  Eq.  (13)  are  the  (absolute  value  of)  difference  and 

sum  of  ko  and  ki  .  jspectively.  The  cross  sections  are  averaged  over  the  initial  magnetic 
substrates  m^  as  indicated  in  Eq.  (14).  The  total  ionization  cross  section  is 

(d^j 

where  E  is  the  incident-electron  energy  and 

E^  =  E-\I\-kll2.  (16) 

The  upper  limit  in  Eq.  (15)  is  set  to  ’Etmnrf^  instead  of  Eimoxj  since  a  close  examination 
[20-22]  shows  that  g(fc2,fci)  and  f(fci,fc2)  describe  the  same  physical  occurrence.  To  avoid  this 
redundancy,  the  upper  limit  in  Eq.  (15)  is  set  to  Emox/2. 

Comparison  of  the  BO  cross  sections  for  ionization  of  the  H  atom  with  those  obtained 
by  the  Bom  approximation  shows  that  inclusion  of  exchange  does  not  greatly  alter  the  cross 
sections  except  at  low  energies.  For  example,  by  reading  the  figure  in  Ref.  18,  we  see  that 
the  ratios  of  the  BO  cross  sections  to  the  Bora  cross  sections  are  about  0.65,  0.75,  0.83,  and 
0.89  at  25,  50,  100,  and  200  eV  respectively. 

B.  Method  of  Distorted  Waves 

In  the  Born  approximation,  plane  waves  e^’^  zmd  e*^*  '  are  used  for  the  incident  and 
scattered  waves.  Because  of  the  interaction  of  the  projectile  electron  with  the  target  atomic 
electrons  and  nucleus,  dep2irture  from  plane  wave  is  expected.  The  method  of  distorted 
wave  (DW)  attempts  to  correct  for  this  departure,  and  has  been  applied  to  electron-impact 
excitations  and  ionization  of  atoms  [23-25].  To  simplify  the  formulation  we  may,  without 
loss  of  generality,  assume  the  direction  of  the  incident  electron  to  be  the  z-direction.  With 
allowance  for  distortion  by  the  target,  the  incident-electron  function  is  expanded  as 

•HKf)  =  (17) 
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while  the  scattered- wave  function  is  expanded  as 


=  41,  (18) 

where  T~^Via{t)  denotes  the  distorted  partial  waves  which  are  the  solutions  of  the 
Schroedinger  equation  for  the  incident  (or  scattered)  electron  in  the  field  of  the  hydro¬ 
gen  atom.  The  precise  form  of  the  differential  equation  for  the  distorted  peirtial  waves  wiU 
be  given  in  Section  2C  when  we  discuss  the  electron-impact  ionization  of  the  oxygen  atom 
[Eq.(42)].  To  obtain  the  DW  scattering  amplitude  we  substitute  Blqs.(l7)  and  (18)  for  e^*" 
and  e’*'  ’''  respectively  in  Eq.  (5),  along  with 


The  DW  direct  amplitude  becomes 

fikxX)  =  -f  1)'/=* 

X  E,c"(fi0,/2m,)c"(fm,/„m„)  Yr^{k)Yi:^,{k) 

X  ^f)»  (20) 

where,  as  defined  by  Condon  and  Shortley  [26], 

c^ihmuhms)  =  [47r/(2A  +  J Yx,m4-m3{r)Yi,mA^)dr,  (21) 

and 

R\kJi,nJ„,kil2,nl)  -  7^fc,Ji(’’)  "Pk^Mdr 

X  j  Pn,u{^)Pni{x)x^dx +r^  J  P„,i.(i)P„i(i)i“'^“'di  .  (22) 

Upon  averaging  over  the  initial  nio  sublevels  and  integrating  over  the  directions  of  the  scat¬ 
tered  electron  and  the  ejected  electron,  we  have  the  DW  cross  sections  (neglecting  exchange) 
as 

(S)  "  fcl*2[(2/o+l)fco]"‘Sm./|/(fcl,  fc2)l*dfeldfc2 

=  16l?kik2[{2l„  +  l)fco]~^£jTniilam3(2fl  -f  1) 

X  Ex[c^(/i0,  l2m2)c\lm,  lom„  =  171  —  7712) 

X  R^{koli,n„lo,kil24nl)]’.  (23) 
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The  exchange  amplitude  can  be  obtained  from  Eq.(6)  by  the  similar  substitutions  for  thi 
incident  and  scatter*  <1  vaves.  However,  such  an  analog  of  the  Bom-Oppenheimer  exchange 
amplitude  leads  to  unreasonably  large  cross  sections.  Therefore,  we  shall  again  resort  to 
Ochkur’s  results  [18],  and  assume  that  f  and  g  are  related  by  Eq.(7).  However,  since  Eq.  (7) 
was  based  on  the  plane- wave  approximation,  the  use  of  this  equation  to  determine  g  from  f  in 
a  DW-type  calculation  must  be  viewed  as  a  further  step  of  approximation.  Nevertheless  tlie 
contribution  to  the  ionization  cross  sections  from  g  is  usually  much  smaller  than  from  f.  Thus 
a  DW  calculation  along  with  an  Ochkur-like  exchange  should  be  a  significant  improvement 
over  the  Bom-Oppenheimer  approximation  for  ionization  calculations.  Following  Eq.  (7), 
we  write 

gif  =  1^0  -  fell V  (*o  - 

=  [kl  +  kl-2Kk,cose,)l[kl~)4) 

=  (a  —  bcosOi),  (2''.) 


where  a  and  b  (^o  +  f^)l{kl  -  fc|)  and  2k„ki/{kl  —  fc|)  respectively,  and  Oi  is  the  angle 
fci  mtikes  with  z-axis  {ko  being  in  the  z-direction). 

Corresponding  to  the  plane-wave  Born-Ochkur  cross  section  we  express  the 

cross  sections  based  on  the  exchange  distorted  wave  method  with  the  Ochkur  approximation 
(EDWO)  as 


kik2  [(2/„  -f- 1)  J  ]f  +  f{a-  bcosO,)  \^dkdk2 

/o  +  A  +  Ij, 


(25) 


where  lo  is  [(a-|-l)*-fb^/3]  times  the  {dalde)DW  of  Eq.  (23),  and 

Ij  =  Oj  ElmX 

X  (ijO,  /2Tn2)  c^  (/m,  lom„  =  m  —  m2) 
X  R  (^^1 ,  ^1^2,  ^^) 


7 


X  c^{l[0,l2m2)c^{lm,ljn„  =  m  —  Tn2) 
X  {k„l[,nolo,kil2,nl) 

X  cP  {l2m2,  ^21712) ,  j  -1,2 


ai  =  — 26(a  +  l),  (27) 

02  =  26V3.  (28) 

Equations  (25)  and  (26)  represent  an  exact  transcription  of  Ochkur’s  formula  to  the  DW 
calculation  of  ionization.  However,  a  much  simpler  expression  is  obtained  in  Eq.(24)  if  we 
make  use  of  the  fact  that  for  excitation  to  an  opticedly  allowed  state  as  well  as  ionization, 
the  scattering  is  strongly  peaked  in  the  forward  direction.  Thus  we  replace  |  by 

(fco  —  kiY  in  Eq.(24),  so  that 


7  =  (30) 

and 

+  l)A^]-‘  (1  f\f  ikuk2)  fdkdk.  (31) 

We  have  made  some  test  cedculations  using  as  example  the  bound-state  excitations  from 
0(3p^  ®P)  to  0[2p^('‘S)3s  ^S]  and  to  0[2p^(^S)3d  ^D].  The  use  of  Eq.(31)  gives  cross  sections 
that  are  larger  than  the  ones  from  Eqs.(25)  and  (26)  by  about  10%  near  the  threshold. 
Beyond  the  threshold,  the  difference  quickly  decreases  to  the  3-5%  ramge  with  increasing 
incident  energy.  In  view  of  such  smcdl  differences,  we  use  this  simplified  version  of  Ochkur’s 
approximation  in  our  calculations. 

C.  Ionization  of  Oxygen  Atoms 

In  this  paper  we  consider  the  ionization  processes; 

e-  -h  O  (2p'‘  ®P)  e-  -h  e-  +  0+  (2p®  "5,  ^D,  ^P)  ,  (32) 

e-  -f  0[(2p^"5)n„/„X]-^e--he-+0+(2p^‘‘5).  (33) 


We  first  determine  the  relevant  orbitals  of  the  oxygen  atom  (in  the  absence  of  the  colliding 
electron)  corresponding  to  the  0(2p^  ^P)  and  0[(2p^  ^S)  ^Loj  bound  states  and  the 
[0'*'(2p®  ^S,*D,*P)  +  e“]  continuum  states  by  the  Hartree-Fock  method  with  a  frozen-core 
approximation  in  which  the  ls,2s,2p  orbitals  are  treated  as  fixed  and  identical  to  the  corre¬ 
sponding  orbitals  for  the  (2p^)^P  state  of  thi  oxygen  atom  as  given  by  dementi  and  Roctti 
[19].  Detailed  procedures  for  calculting  the  no^o-orbitals  and  the  continuum-state  orbitals 
have  been  described  previously  [27,28].  Analogous  to  the  notation  in  Section  2A,  we  use 
<^o(^  for  the  active  electron  (2p  or  no^-o)  of  the  initial  target  sta.te,  and  4>j{k2ir)  for  the 
ejected  electron.  The  latter  is  a  superposition  of  anguleir  momentum  states  of  different  i,n\ 
values  and  behave  at  large  distance  like  a  pzirticle  of  momentum  hk  as  indicated  in  Ref.  28. 
It  is  worth  noting  that  the  orthogonality  condition  is  incorporated  into  the  differential  equa¬ 
tion  for  the  continuum-state  orbitals  [28],  thereby  obviating  the  arbitrary  orthogonalization 
procedure  at  a  later  stage  of  the  calculation. 

Ionization  is  produced  by  the  interaction  of  the  colliding  electron  with  the  0  atom 

^in.  =  --  +  S’  |ri-r^|-S  (34) 

ri 

where  the  nuclear  charge  Z  is  8,  and  we  labeled  the  colliding  electron  as  1.  Since 
is  independent  of  spin,  we  construct  the  spin  eigenfunctions  of  the  colhding  system.  For 
example,  in  the  case  of  Eq.  (32)  we  have  for  the  initial  state 

^o(5, M|xi,i2,— ^b)  =  .AIl^C(l/2,  p,  1,  M-p;  SM) 

<f)c  (1/2,  pjxx)  $<>  (1,  M  -  p|i2, ...,  19) ,  (35) 

where  C(ji,  mi,  j2,  m2;  jm)  is  the  Clebsch-Gordan  coefficient,  denotes  the  wave  functions 
of  the  oxygen  atom  in  the  (2p^)®P  state,  and  the  colliding  electron  is  represented  by  (f)^  which 
will  be  described  later.  We  use  Xj  to  represent  the  spatial-spin  coordinates,  and  A  for  the 
antisymmetrization  operator.  To  construct  the  final-state  wave  function  for  Eq.(32),  we  first 
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couple  the  O'*' -ion  function  (#''■)  with  the  spin-orbital  of  the  ejected  electron  to  form 
the  intermediate  states,  i.e., 


{1/2,  (i,  S,M,) 

4>ej  (l/2,  fcajica)  (3/2,  Afj  -  /i|i3, 19) ,  (36) 


which  are  then  coupled  with  the  colliding  electron  function  <f>c  to  form  'I'/(S'M'  1 
The  scattering  amplitudes  are  obtained  by 

=  (l/27r)  J  ^){SM\xu...,£9)Hiru^o{SM\  xr,...,x^)dx^...dx^,  (37) 

which  axe  diagonal  in  (SM)  and  independent  of  M.  As  a  consequence  of  and  $  /  being 
antisymmetrized,  includes  both  the  direct  f  and  exchange  g  amplitudes. 

Next  we  sum  the  cross  sections  over  S  and  M,  and  divide  by  2(2Si-l-l),  the  number  of 
the  possible  initial  spin  states.  Si  being  the  spin  of  the  initial  atomic  state,  so  that  the 
results  correspond  to  ionization  by  a  beam  of  unpolarized  electrons.  By  a  straight-forward 
calculation,  we  find 

J  =  (A;iV*o)£5[(25-hl)/2(2S, +l)]/|.F^|Vfcidfc2  (38) 

=  {k^k^/K)  I  (1/1^  -fg+  \g\^)  dkdk,  (39) 

/  =  -2K-^CFpJ<f>o{r^e^^U}{i^,r)df,  (40) 


where  Cfp  stands  for  the  coefficients  of  fractional  parentage  that  arise  due  to  the  equivalent 
2p  orbitals  in  process  (32),  and  is  numerically  equal  to  ^4/3,  ^5/3  ,and  1  respectively  for 
the  '‘S,  *D  and  final  states  of  0"*'[29].  For  process  (33)  the  issue  of  fractional  parentage 
does  not  £inse,  i.e.,  Cfp=1.  The  g  amplitude  is  related  to  f  through  Eq.(7)  or  Eq.(29). 
Notice  that  there  are  no  superscipts  on  the  left-hand  side  of  Eq.(38),  since  the  cross  section 
corresponds  to  ionization  by  unpolarized  electrons. 
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To  calculate  the  BO  cross  sections,  4>e{f)  is  taken  as  the  incident  plane  wave  exp{iko-'i^ 
in  and  as  exp(i^i-f)  in  'if,  and  Eq.(7)  is  used  to  express  g  in  terms  of  f.  This  gives 

'  '  BO 

X  \  <nX\K\£>\^K-^dK.  (41) 

For  the  EDWO  calculations,  <f)c  takes  on  the  more  complicated  form  as  given  by  Eqs.(17) 
and  (18).  In  these  equations  the  distorted  partial  waves  are  the  solution  of 

-V(r)  +  A:^]  Vu  (r)  =  (r)  (r)  +  (r) ,  (42) 

where  the  direct  potential  V(r)  is  due  to  all  eight  atomic  electrons  in  the  oxygen  atom  foi 
the  appropriate  initizd  state  in  process  (32)  or  (33).  Similarly,  represents  the 

exchange  interaction  with  all  the  target  electrons.  The  last  term  in  Eq.(42)  is  to  ensure  that 
V  is  orthogonal  to  all  the  bound  ni  orbitals  via  the  Lagrange  multiplier  Cn»j.  The  dettuls 
are  further  described  in  Ref.  28.  The  distorted-wave  functions  of  the  projectile  electron  have 
the  asymtotic  form  (as  r  oo)of 

Vu{,r)  ~  k~^sin  ^kr  —  ^/tt  -f-  ,  (43) 

where  »s  the  phase  shift  due  to  the  potentials.  As  explained  earlier,  for  the  distorted- wave 
calcultion  we  use  the  simplified  version  of  the  Ochkur  exchange  as  given  in  Eq.(29).  The 
EDWO  cross  sections  then  become 

l&7rkik2CFp[ko  {2lo  +  (2/i  +  1)  [1  -  7  -H  7^] 

(fiO,  127712)  {lm,lomo  =  m  —  71x2) 

R^{Kli,nM,k^l2,nl)]\  (44) 

where  7  is  defined  in  Eq.(30).  The  total  cross  sections  (Tbo  and  (Tedwo  are  then  obtained 
by  integrating  over  e  as  indicated  in  Elq.  (15). 
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3.  RESULTS 

In  order  to  check  the  numerical  procedures  we  have  computed  the  ionization  cross  sections 
of  H(2s)  using  the  numerically  tabulated  functions  for  both  the  bound  2s  state  and  the 
continuum  (Coulomb)  functions.  In  this  test  calculation  we  included  the  angular  momenta 
of  the  ejected  electron  ^  =  0  -  8.  We  were  able  to  reproduce  the  cross  sections  computed  by 
Prasad  [30]  within  4%  for  incident-electron  energy  from  threshold  to  200  eV. 

In  general  a  greater  number  of  angulu  momenta  I  of  the  ejected  electron  axe  needed 
to  achieve  convergence  as  the  incident-electron  energy  is  increased.  Also,  the  higher  the 
quantum  state  from  which  ionization  takes  place,  the  more  Vs  are  needed.  In  this  paper 
we  used  ^  =  0  to  /  =  30  in  the  BO  calculations  for  ionization  of  the  oxygen  atoms  in  the 
ground  state  and  excited  states  (noio)  up  to  no  =  6,  and  used  ^  =  0  to  f  =  10  for  the  ED  WO 
calculations  which  cover  the  groimd  state  and  excited  states  with  no  =  3. 

The  integration  with  respect  to  K  in  Eq.  (13)  is  carried  out  numerically.  The  smallest 
step  size  AK  is  4  x  10~'*,  and  the  step  size  is  doubled  after  10  quadrature  points,  except  for 
the  eleventh  (last)  region  in  which  there  are  30  points,  for  a  total  of  130  points.  The  rather 
small  AK  in  the  beginning  is  to  take  care  of  the  steeply  rising  differential  cross  sections  in 
the  forward  direction  (i.e.,  toward  small  K)  of  the  dipole-allowed  tr^msitions,  to  which  all 
ionization  processes  belong.  Sinularly,  with  regard  to  e-integration  in  Eq.  (15),  we  used  the 
initial  Ae  =  7.8125  xl0~^  a.u.,  and  doubled  it  after  each  region.  There  are  eight  regions 
with  10,  10,  10,  10,  40,  20,  20,  and  80  quadrature  points  each. 

A.  Ionization  of  0(2p'*  ^P)  to  Form  0'^(2p®  ^S,  *D,  *P) 

We  have  computed  the  Bom  (B)  and  BO  cross  sections,  in  the  remge  of  17.5  -  500  eV  of 
the  incident  electron  energy,  of  ionizing  the  ground  state  0(2p^  ®P)  atoms  and  leaving  the 
0+(2p^)  ions  in  the  three  states,  ^S,  ’D,  and  *P.  The  sum  of  these  three  cross  sections,  which 
corresponds  to  ionization  of  the  ground-state  oxygen  atoms  to  form  0‘''(2p^)  ions,  is  shown 
in  Table  1.  We  have  also  applied  the  EDWO  method  as  well  as  the  method  of  distorted 
waves  without  exchange  (DW)  to  calculate  the  cross  sections  at  incident  energies  1 7.5  -  200 
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eV  and  the  results  aie  included  in  Table  1.  Comparison  of  the  B  with  BO  cross  sections  and 
of  the  DW  with  ED  WO  cross  sections  shows  that  the  effect  of  electron  exchange  is  to  reduce 
the  cross  sections.  At  25  eV  this  reduction  amounts  to  11%  for  B  and  18%  for  DW,  whereas 
at  200  eV,  the  BO  is  only  5%  below  the  B  cross  section,  jmd  the  difference  between  EDWO 
and  DW  is  even  smaller.  Comparison  of  the  BO  and  EDWO  columns  in  Table  1  underscores 
the  importance  of  the  distortion  of  the  projectile  at  low  energies.  The  EDWO  cross  sections 
are  smaller  thain  the  BO  cross  sections  at  25  and  50  eV,  but  the  order  is  reversed  at  75-200 
eV.  Furthermore  the  EDWO  and  BO  cross  sections  tend  to  merge  at  increasing  energies. 

The  ionization  cross  sections  of  0(^P)  atom  have  been  calculated  by  several  investigators. 
Peach  [11,12]  computed  ionization  cross  sections  resulting  in  0‘‘'('*S,  ^D,  ^P)  from  the  ground 
state.  In  the  work  of  Peach  [11]  the  function  of  Ref.  17  was  used  for  the  2p  orbital,  and 
undistorted  Coulomb  functions  of  unit  nuclear  charge  (hydrogenic  Coulomb  functions)  were 
used  for  the  ejected  electron.  The  Bom  and  Bom-Ochkur  methods  were  used,  thus  Peach’s 
work  is  similar  to  the  present  one  except  for  the  use  of  the  hydrogenic  Coulomb  functions. 
Later  the  author  [12]  discovered  an  error.  As  a  result  the  cross  sections  [11]  are  to  be 
reduced  by  a  factor  of  two.  McGuire  [13]  computed  the  Bom  cross  sections  by  using  the 
wave  functions  based  on  the  Herman- Skillman  potential.  Omidvar  et  al.  [14]  used  the 
scaled  hydrogenic  functions  to  compute  the  Bom  cross  sections.  Kazaks  et  al.  [15]  used 
the  potential  of  the  independent-particle  model  [31]  to  determine  the  wave  functions  from 
which  the  Bom  cross  sections  were  obtained.  Burnett  and  Rountree  [16]  reported  a  Born- 
approximation  calculation  using  a  set  of  very  refined  target  wave  functions.  The  latter  four 
works  [13-16]  do  not  include  the  effect  of  electron  exchange.  In  Table  2  we  compare  our  results 
with  those  of  the  previous  works.  The  columns  labeled  as  BAR  (Burnett  and  Rountree  [16]), 
OKS  (Omidvar  et  al.  [14]),  and  KGG  (Kazaks  et  al.  [15])  are  our  best  reading  of  the  figures 
in  the  respective  papers.  The  column  PEA  (Peach  [11])  shows  the  interpolation  of  Table  4  of 
Ref.  11  divided  by  2  as  corrected  by  Ref.  12.  Likewise,  the  column  MCG  (McGuire  [13])  is  an 
interpolation  of  Table  4  of  Ref.  13.  The  target  wave  functions  used  by  Burnett  and  Rountree 


13 


TABLE  1.  Electron-impact  ionization  cross  sections  (in  10“^^  cra^)  of 
ground-state  oxygen  atoms  resulting  in  O'*"  ions  in  the  (2p^)^S,  ^D,  and  states 
calculated  in  the  present  work  by  using  the  method  of  exchange  distorted  waves 
with  Ochkur*s  approximation  (EDWO);  by  the  method  of  distorted  waves  with  no 
exchange  (DW);  by  the  Born-Ochkur  (BO)  approximations  and  by  the  Born  (B) 
approximation  with  no  exchange. 


E(eV) 

EDWO 

DW 

BO 

B 

25 

0.171 

0.210 

0.326 

0.367 

50 

0.952 

1.027 

0.996 

1.123 

75 

1.287 

1.347 

1.252 

1.367 

100 

1.360 

1.405 

1.293 

1.411 

125 

1.337 

1.369 

1.272 

1.274 

150 

1.276 

1.302 

1.223 

1.310 

200 


1.134 


1.149 


1.091 


1.152 


TABLE  2.  Electron-impact  ionization  cross  sections  (in  10"^^  cm^)  of  oxyget. 
atoms  to  form  O'*'  ions  in  the  (2p^)^S,  -D,  and  states  calculated  by  the  method 
of  exchange  distorted  waves  with  Ochkur's  approximation  (EDWO)  and  by  the 
Born-Ochkur  approximation  (BO)  in  this  work  at  different  electron  energies  E  (in 
eV).  Included  for  comparison  are  the  cross  secstions  from  the  earlier  worki.-  by 
Burnett  and  Rountree  (BAR)  [16],  by  Peach  (PEA)  [11,12],  by  McGuire  (MCG),  [131 
by  Omidvar,  Kyle,  and  Sullivan  (OKS)  [14],  and  by  Kazaks,  Ganas,  and  Green  (KGSj 
[15]  as  explained  in  the  text.  The  last  five  sets  of  cross  sections  except  PEA 
were  calculated  by  means  of  the  Born  approximation  (B)  rather  than  the 
Born-Ochkur  approximation  (BO). 

This  work 


E(eV) 

EDWO 

BO 

BAR(B) 

PE A (80) 

MCG(B) 

OKS(B) 

KGG(B) 

25 

0.171 

0.326 

0.38 

0.192 

0.55 

0.55 

50 

0.952 

0.996 

1.12 

0.625 

1.15 

1.42 

2.2 

75 

1.287 

1.367 

1.36 

0.867 

1.35 

1.77 

2.2 

100 

1.360 

1.293 

1.38 

0.949 

1.48 

1.85 

2.1 

125 

1.337 

1.272 

1.30 

0.989 

• 

CO 

1.76 

1.9 

150 

1.276 

1.223 

1.26 

0.991 

1.44 

1.65 

1.8 

200 

1.134 

1.091 

1.14 

0.927 

1.26 

1.51 

1.5 

300 

0.860 

0.790 

0.97 

1.20 

1.2 
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are  more  refined  than  ours.  For  the  initial  target  state,  they  used  a  six-term  Cl  wave  function 
as  opposed  to  our  one-configuration  Hartree-Fock  work.  To  find  the  continuum  functions 
for  the  ejected  electron,  Burnett  and  Rountree  solved  the  close-coupling  equations  iciuding 
the  ^S,  *D,  *P  term  manifold  of  the  0'''(2p®)  configuration  whereas  in  our  case  the  ion 
is  taken  to  be  in  the  (2p*)^S  state.  Comparing  the  cross  sections  of  Burnett  and  Rountree 
(BAR),  which  were  calculated  by  means  of  the  Bom  approximation,  listed  in  Table  2  with 
our  cross  sections  calculated  by  the  same  method  given  in  the  last  column  of  Table  1 ,  we 
see  an  agreement  within  4%.  This  supports  the  accuracy  of  our  wave  functions.  The  cross 
sections  listed  imder  PEA  (BO)  in  Table  2,  which  were  determined  by  Peach  using  the  BO 
method,  are  substantially  smaller  than  our  BO  results  (about  30%).  We  believe  that  this 
is  largely  due  to  the  different  continuum  wave  functions  (for  the  final  target  states)  used  to 
obtain  the  two  sets  of  cross  sections  since  the  initial  target  state  wave  functions  are  very 
similar  in  Peach’s  and  our  calculations.  In  the  cases  of  the  MCG,  OKS,  and  KGG  calculations 
(Table  2),  the  wave  functions  for  both  the  initial  and  final  target  states  are  different  from 
the  ones  used  in  the  present  work  and  their  cross  sections  differ  from  our  Bom  cross  sections 
to  varying  extents. 

Measurements  of  the  ionization  cross  sections  of  0  atom  were  reported  by  Fite  and 
Brackmann  (FB)  [4].  In  a  cross-beam  experiment,  they  determined  the  ratio  of  the  cross 
sections  for  ionizing  atomic  oxygen  and  molecular  oxygen,  o’(0''')/<r(02‘*').  Then  utilizing 
the  molecular  ionization  data  of  Tate  and  Smith  [8],  they  obtained  <r(0‘^).  In  a  similar 
experiment  Rothe  et  al.[5]  measured  ionization  cross  sections,  and  their  results  are  in  good 
zigreement  with  those  of  FB  [4].  In  the  experiment  of  Brook,  Harrison,  and  Smith  (BHS) 
[6],  the  charge-exchange  process  was  used  as  the  sotirce  of  atomic  beam  so  that  the  absolute 
cross  sections  were  determined  directly.  Zipf  [7],  using  the  new  data  of  02'''  ionization  of 
Mark  [9],  revised  the  cross  sections  of  FB.  The  “revised”  cross  sections  are  smaller  than  the 
uncorrected  ones  by  a  substantial  margin  below  100  eV,  but  above  100  eV  the  difference 
is  only  about  10%.  Zipf  [7]  also  determined  the  cross  sections  o’(0‘*')  by  normalizing  his 
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measurements  to  the  data  of  BHS  and  the  revised  data  of  FB  in  the  100-300  eV  range.  The 
shape  of  the  ionization  function  (a  plot  of  cross  sections  vs  incident  energy)  of  Zipf  [7j  agrees 
quite  well  BHS  and  also  with  FB  above  100  eV. 

In  Figure  1  we  compare  the  BO  and  EDWO  cross  sections  of  this  work  with  the  exper 
imental  data  of  BHS  and  of  FB  as  revised  by  Zipf.  [7].  The  data  of  Zipf  are  not  included 
because  they  are  quite  close  to  BHS.  The  agreement  between  the  present  calculation  and 
experiments  is  seen  to  be  quite  good,  especially  at  energies  above  70  eV.  In  the  70- 1.')!)  eV 
rzmge  the  difference  is  only  about  3%,  and  at  200  eV  the  present  Ccdculation  is  about  10% 
smaller  thzm  the  experimental  value.  A  greater  difference  is  found  at  lower  energies;  for 
example,  at  50  eV  our  value  is  about  20%  smaller  than  the  experimentad  value.  Below  70 
eV  there  exists  a  substantial  difference  between  the  two  sets  of  experimental  data  shown  in 
Figure  1.  Since  the  BHS  data  came  from  a  more  recent  experiment,  we  use  them  as  the 
primary  source  of  experimental  data  for  comparison  with  our  calculations.  In  comparing 
our  calculation  with  experiments,  one  must  remember  that  the  experimental  cross  sections 
include  the  effects  of  autoionization  as  well  as  the  direct  ionization,  whereas  only  the  lat¬ 
ter  is  considered  in  the  theory.  While  Brook  et  al.  [6]  indicated  that  their  experimental 
data  show  no  discernible  evidence  of  autoionization,  detailed  examinations  of  the  ionization 
yield  and  optical  emission  of  the  autoionizing  states  [32-34]  show  that  the  influence  of  au¬ 
toionization  on  the  measured  cross  sections  must  be  taken  into  consideration.  In  particular 
autoionization  via  the  O  [2s*2p^(*P)38  "  ®P“]  and  O  [2s2p®  ^P°]  states  are  foimd  to  contribute 
significantly  to  the  observed  ionization  cross  sections.  Since  the  ^P®  symmetry  cannot  be 
realized  from  the  [2p®  ^S]  ion  plus  an  sulditional  (unbound)  electron,  the  autoionization 
rates  of  the  above-mentioned  ®P®  states  are  much  reduced  and  become  comparable  to  the 
radiative  emission  rates.  Dehmer  et  al.  [32]  experimentally  determined  the  ratio  of  the 
autoionization  rate  to  the  emission  rate  to  be  nearly  equal  to  unity  for  both  the  3s"  ^P° 
and  2s2p®  ®P®  states,  i.e.,  half  of  the  population  in  these  states  decay  by  autoionization 
and  the  other  half  by  emission.  Zipf  and  Kao  [33],  from  their  experimental  emission  cross 
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sections  and  the  branching  ratios  of  Dehmer  et  al.  [32],  conclude  that  autoionization  con 
tributes  substantially  to  the  production  of  0"^(^S)  ions  by  electron  impact.  In  a  more  recent 
electron-impact  energy  loss  experiment,  Vaughan  and  Doering  [34]  determined  the  excita¬ 
tion  cross  sections  for  several  states  of  atomic  oxygen  including  the  two  autoionizing  ^P*’ 
states  mentioned  earlier  and  the  2p^(^D)4d'  state  at  incident-electron  energies  betwei-n 
30  and  200  eV.  The  excitation  functions  for  these  states  are  seen  to  peak  at  50  eV  [31’. 
Applying  the  branching  ratios  of  Dehmer  et  al.  referred  to  ecirlier  to  the  excitation  cross 
sections  of  Vaughan  and  Doering  for  the  3s"  £ind  2s2p^  ^P“  states,  we  can  determine  the 
contribution  to  the  ionization  cross  sections  due  to  autoionization  through  the  3s"  ^P"  and 
2s2p®  ^P°  states.  The  excitation  cross  section  for  the  2p^(^D)4d'  ^P“  state  is  mucli  smaller 
than  those  of  the  other  two  ^P"  states  mentioned  above;  thus  even  if  we  assume  that  the 
2p®(^D)4d'  ^P"  state  decays  by  autoionization  entirely,  this  autoionization  channel  mcikes 
a  much  smaller  contribution  to  the  measured  ionization  cross  section  compared  to  the  3s" 
^P®  and  2s2p®  ^P®  channels.  Indeed  Dehmer  et  al.  [32]  indicated  that  the  2s*2p^(^D)3d'’ 
^P"  state  decays  primarily  by  autoionization.  Thus,  we  assume  a  100%  autoionization  decay 
for  the  2s^2p^(^D)4d'  ^P“  state  in  order  to  determine  the  contribution  from  this  channel  to 
the  ionization  cross  section.  From  this  analysis  we  obtain  the  autoionization  contributions 
to  the  observed  ionization  cross  sections  as  0.068,  0.151,  0.090,  0.071,  and  0.048  x  10"'*® 
cm^  at  30,  50,  100,  150,  and  200  eV  respectively.  When  this  “correction”  is  made  to  the 
measured  ionization  cross  sections  of  BHS,  the  agreement  between  the  present  theory  and 
experiment  is  much  improved  at  lower  energies  as  shown  in  Figure  2.  For  instance  at  50  eV 
our  theoretical  vzdue  is  only  5%  smaller  than  the  experimental  data,  whereas  at  100  eV  our 
value  is  now  8%  smaller  and  at  200  eV  6%  larger  than  the  experimental  counterparts.  At  32 
eV,  the  lowest  energy  for  which  we  can  make  the  autoionization  correction,  our  theoretical 
value  is  36%  lower  them  the  experimental  cross  section. 
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B.  Ionization  of  0[2p^(^S)noA>  ^Lo]  to  Form  0'''(^S) 

We  have  calculated  the  cross  sections  for  ionizing  the  excited  oxygen  .items  0[2p  *(  ’S)ii,/„ 
^Lo]  by  electron  impact.  In  Figures  3  and  4,  we  show  the  BO  eind  EDWO  cross  sections  for 
n^^o  =  3s,  3p,  zmd  3d.  In  each  case  the  pealc  of  the  ionization  function  occurs  immediately 
above  the  threshold  (8.0  eV  for  3s,  6.0  eV  for  3p,  and  4.0  eV  for  3d),  in  marked  contrast  to 
ionization  of  the  ground  state  where  the  peak  cross  section  occurs  at  about  100  eV.  Such  a 
characteristic  distinction  between  ionization  from  the  ground  state  and  from  excited  states 
is  also  seen  for  the  H  and  He  atoms.  The  ionization  functions  for  H(2s)  and  H(2p)  p<!ak 
around  8-13  eV  [30],  whereas  the  maximum  of  the  H(ls)  ionization  function  occurs  within 
80-100  eV.  Similarly  the  ionization  functions  for  He(2^S)  and  He(2^S)  peak  around  10  eV 
[35],  as  opposed  to  the  peak  position  in  the  80-100  eV  range  for  He(ls)  [36].  As  in  the  case 
of  ionization  from  the  groxmd  state,  the  EDWO  cross  sections  are  substantially  smedler  than 
the  BO  cross  sections  at  low  energies.  However,  the  two  sets  of  cross  sections  become  much 
closer  to  each  other  with  increasing  energy  as  can  be  seen  in  Table  3. 

We  have  extended  the  BO  calculation  to  n^Aj  up  to  (6,5).  The  general  shape  of  the 
ionization  function  remains  the  same  except  for  a  small  inward  shift  in  the  peak.  We  have 
fitted  the  ionization  function  for  each  Ugio  from  7.5  eV  to  500  eV  to  the  relation 

<T  =  c\E~^lnE  -h  C2E~^  -f  C3E~^,  (45) 

where  the  incident-electron  energy  E  is  in  eV  and  the  cross  section  in  cm*.  The  numerical 
values  of  Ci,  C2,  and  C3  are  given  in  Table  4. 

Finally  it  should  be  recalled  that  in  the  present  calculation  the  initial  states  are  triplet 
states.  Although  we  did  not  make  a  separate  set  of  cedculations  for  the  case  of  quintet 
initial  states,  we  expect  the  latter  cross  sections  to  be  rather  close  to  the  former  set.  In  the 
radiative-recombination  calculation  [28]  involving  the  above  two  sets  of  triplet  and  quintet 
states,  the  difference  in  cross  sections  were  found  to  be  about  15%  for  n,,  =  3,  10%  for  n^  = 
4,  and  less  still  for  higher  n^,. 
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TABLE  3. 


Comparison  of  the  BO  and  EDWO  cross  sections  for  ionization  of  oxygen 
atoms  in  the  [2p^(**S)nQj^]3LQ  states.  For  each  incident-electron  energy  E,  there 
are  two  entries:  the  first  is  the  BO  and  the  second  the  EDWO  cross  sections  (in 
10"^^  cm^). 


E(eV) 

3s 

3p 

3d 

2.5 

7.49 

3.90 

5.0 

0.659 

2.80 

8.62 

0.189 

1.S8 

7.36 

7.5 

1.19 

2.75 

7.14 

0.700 

2.21 

6.37 

10.0 

1.14 

2.43 

5.93 

0.818 

2.08 

5.34 

12.5 

1.02 

2.12 

5.03 

0.815 

1.88 

4.54 

15.0 

0.917 

1.88 

4.35 

0.774 

1.69 

3.93 

17.5 

0.827 

1.68 

3.B5 

0.724 

1.53 

3.49 

20.0 

0.752 

1.51 

3.44 

0.674 

1.39 

3.00 
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TABLE  4.  Values  of  the  coefficients  in  the  Bethe-type  expansion  [Eq,  (45)]  of 

the  ionization  cross  sections  (in  cm^)  for  the  oxygen  atoms  in  the 

[2p^( states  as  a  function  of  the  incident-electron  energies  (in  eV) 


"o-^o 

=1 

=2 

^3 

3.0 

1.135[-15]^ 

1.458C-14] 

-5.797[-14] 

4,0 

2.582[-15] 

4.172C-14] 

-7.230C-14] 

5,0 

3.760[-15] 

8.517L-14] 

-8.563E-14] 

6,0 

1.473C-15] 

1.480C-13] 

-1.096E-13] 

3.1 

1.54U-15] 

3,044[-14] 

-9.684C-14] 

4,1 

7.842[-15] 

5.170C-14] 

-5.615E-14] 

5.1 

1 .522[-14] 

9.016C-14] 

-5.832E-14] 

5.1 

1.709[-14] 

1.458C-13] 

-5. 230 [-14] 

3,2 

4.265C-15] 

6.280C-14] 

-13.380C-14] 

4.2 

1.354C-14] 

8.649C-14] 

-6.529E-14] 

5.2 

2.437C-14] 

1.202C-13] 

-3.767E-14] 

6,2 

1.627C-14] 

2.028C-13] 

-6.385C-14] 

4,3 

7.23U-15] 

1,084[-13] 

-9.185C-14] 

5.3 

1.688C-14] 

1.444C-13] 

-6.222C-14] 

6,3 

1.79U-14] 

2.013[-13] 

-5.938C-14] 

5,4 

7.499C-15] 

1.696[-13] 

-8.738[-14] 

6,4 

1.088C-14] 

2.239C-13] 

-7.515E-14] 

6,5 

7.535C-15] 

2.309C-13] 

-7.217E-14] 

^Number  inside  the  brackets  indicates  the  power  of  10. 
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4.  SUMMARY  AND  CONCLUSIONS 

We  have  calculated  the  electron-impact  ionization  cross  sections  of  producing  the  O'*^ 
ions  (in  the  2p®  configuration)  from  the  ground  state  0(®P)  atom  by  the  BO  and  by  ED  WO 
approximations.  For  cx>mparison  similar  calculations  neglecting  electron  exchange  have  been 
performed  (the  B  and  DW  approximations).  The  continuum  wave  functions  of  the  ionized 
electron  were  computed  by  the  HF  method  with  full  allowance  for  the  Coulomb  and  exchange 
interactions  with  the  remaining  electrons  and  the  nucleus  of  O"*" . 

The  EDWO  gives  significantly  smaller  cross  sections  than  does  BO  at  low  electron  ener¬ 
gies.  For  example,  at  30  eV  the  EWDO  cross  section  is  smaller  by  15%,  and  at  lower  energy 
smaller  still  by  a  greater  percentage.  However,  between  70-200  eV,  the  difference  is  about 
10%  or  less,  and  the  two  sets  of  cross  sections  tend  to  merge  at  high  incident  energy. 

By  accounting  for  the  distortion  of  the  projectile  electron  by  the  target  and  using  the 
proper  continuum  functions,  we  made  a  significant  improvement  over  the  previous  czdcula- 
tions,  which  are  within  the  stage  of  the  B  or  BO  approximation  and,  in  most  czises,  with 
the  hydrogenic  Coulomb  functions,  so  that  we  can  make  a  quantitative  comparison  with 
the  experimental  measurements  [4-7].  In  the  energy  range  of  50-200  eV  the  present  EDWO 
cross  sections  agree  with  the  experimental  measurements  (with  correction  for  autoionization 
contributions)  within  about  10%  or  less.  Below  50  eV  the  discrepancy  becomes  lEUger.  For 
instance  at  32  eV  our  EDWO  cross  section  is  36%  below  the  experimental  value.  One  may 
speculate  that  the  larger  discrepancy  at  low  energies  may  be  partly  due  to  a  possible  increase 
in  the  percentage  uncertainty  of  the  experimental  measurements  on  account  of  the  cross  sec¬ 
tions  being  much  smaller  and  varying  steeply  with  energy.  It  should  also  be  pointed  out  that 
in  our  BO  and  EDWO  cedculations,  the  polarization  of  the  target-state  wave  functions  by 
the  incident  electron  is  neglected.  This  target  polarization  effect  is  more  important  at  very 
low  incident-electron  energies,  and  may  be  partly  responsible  for  the  discrepancy  between 
theory  and  experiment  at  very  low  energies.  We  may  also  add  that  at  energies  above  37 
eV,  ionization  may  take  place  through  channels  in  which  the  O'*"  ion  is  initially  formed  in 
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2U1  excited  configuration  such  as  2p*3s.  However,  contributions  from  these  channels,  whicli 
involve  two  active  electrons  and  are  not  included  in  our  calculations,  to  the  ionization  cross 
sections  are  unimportant  when  we  are  discussing  accuracy  on  the  level  of  10%. 

We  have  computed  BO  cross  sections  of  ionizing  Rydberg  oxygen  atoms  in  the 
[2p^(^S)no^o]^Lo  for  (no^o)  =  (3,0)  through  (rio^o)  =  (6,5).  These  cross  sections  are  sharply 
peaked  within  a  few  eV  above  the  threshold  in  contrast  to  the  ionization  function  of  the 
ground  state,  which  shows  a  broad  peak  around  100  eV.  Similar  characteristic  distinction 
between  ionization  of  the  ground  state  and  excited  states  has  also  been  reported  for  the  H 
2Uid  He  atoms.  For  comparison  we  have  used  the  EDWO  approximation  to  czJculate  the  ion¬ 
ization  cross  sections  for  the  Rydberg  oxygen  atoms  with  no  =  3.  The  EDWO  cross  sections 
are  significantly  smaller  than  their  BO  counterparts  at  low  incident  energies  but  above  15 
eV  or  so  the  two  sets  become  quite  close  to  each  other. 

To  summarize  the  large  amount  of  the  BO  ionization  cross  section  data  for  the  Rydberg 
oxygen  atoms,  we  present  the  cross  sections  in  the  form  of  a  Bethe-type  expansion.  Detailed 
information  about  ionization  processes  involving  excited  oxygen  atoms  is  important  in  studies 
of  radiation  in  the  upper  atmosphere  and  in  modeling  and  diagnostics  of  plasmas. 
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PART  2 


THREE-BODY  RECOMBINATION  OF  ELECTRONS  WITH  ATOMIC 

OXYGEN  IONS 

The  three-body  recombination  of  electrons  with  atomic  oxygen  ion  to  form  an  excited 
oxygen  atom  in  the  2p*(^S)n/  state, 

e~  4-  e”  +  O'*"  [2p^('*5)j  — ♦  e~  -t-  0  ^2p®(^5)nZ  ,  (46) 

is  the  inverse  process  of  electron-impact  ionization, 

e-  +  0  [2p^CS)nl  ^  e-  4-  e'  -h  0+  [2p®(^5)]  (47) 

discussed  in  Part  1.  With  the  electron-impact  ionization  cross  sections  at  various  incident 
energies  obtained  in  Part  1,  we  can  calculate  the  ionization  rate  coefficients  and  the  three- 
body  recombination  constants. 

Consider  an  oxygen  atom  initially  in  the  2p®(^S)nf  ®L  state,  which  we  designate  as  state 
“p”,  colliding  with  an  electron  of  velocity  v  and  resulting  in  ionization  with  the  ejected 
electron  in  a  continuiun  state  “c”.  If  we  denote  the  cross  section  for  this  process  as  crp(v), 
the  ionization  coefficient  at  temperature  T  corresponding  to  the  p— >c  ionization  process  is 

K{p-*c]T)=f  v<Tp{v)f{v,T)dv,  (48) 

^V0 

where  /(v,  T)  is  the  Maxwell-Boltzmaim  distribution  function, 

/(v,  T)  =  Ai7{ml2wkTf'hxp{-mv^l2kT),  (49) 

and  Vo  is  the  threshold  incident  velocity  for  the  p  c  ionization.  Since  we  have  already 
determined  the  ionization  cross  section  a’p(v)  in  Part  1,  the  ionization  rate  coefficient  can  be 
calculated  as  a  function  of  temperature  from  Eq.  (48). 

The  three-body  recombination  constant  is  denoted  by  K(c  — +  p;  T)  and  is  related  to  the 
ionization  coefficient  through  the  equilibrium  constant  K^iiT)  as 

K{c  -  p;  T)  =  K{p  c;  T)IK^{T\  (50) 
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Ka,{T)  =  2{g+lgp]{2-KmkTlhf/’^exp{-Ip/kT),  (51 ) 

where  and  gp  are  the  degeneracies  of  the  atomic  ion  and  the  state  “p”  of  the  atom 
respectively,  and  Ip  is  the  ionization  energy  of  the  state  “p”.  Note  that  the  ionization  rate 
coefficient  is  customarily  expressed  in  units  of  cm^/s.  The  equilibrium  constant  is  in 
(length)"^,  so  that  the  three-body  recombination  rate  constant  K[c  P\T)  is  expressible 
in  units  cm®/s. 

The  ionization  rate  coefficients  of  the  ground-state  oxygen  atom  are  computed  with  the 
cross  sections  obtained  by  the  exchange-distorted-waves  method  with  Ochkur’s  approxima¬ 
tion  (EDWO)  as  described  in  Part  1.  The  results  are  shown  in  Figure  5.  For  the  oxygen 
atoms  in  triplet  excited  states  2p^(^S)n^  (with  n  >  3),  the  ionization  rate  coefficients  are 
czilculated  with  the  Bom-Ochkur  cross  sections  and  the  results  are  shown  in  Figures  6-9. 
Although  we  have  not  computed  the  cross  sections  for  the  quintet  states,  i.e.,  2p^('‘S)nf 
then  cross  sections  are  expected  to  be  very  close  to  the  cross  sections  for  the  triplet  coun¬ 
terparts.  In  the  radiative  combination  process,  the  cross  sections  of  capture  into  the  triplet 
and  quintet  states  are  quite  close  [1]. 

The  three-body  recombination  constants  have  been  calculated  from  Eq.  (50)  and  the 
results  are  shown  in  Figures  10-14.  We  found  that  the  temperature  dependence  for  the 
three-body  recombination  constants  for  the  excited  states  of  oxygen  can  be  fitted  to  the 
form 

K{c-*p-,T)  =  aT-^  (52) 

in  two  regions  of  electron  temperature  T;  5x10^  -  2x10^  ®K  and  2x10*  -  2x10®  "K.  The 
parameters  a  and  ij  are  shown  in  Table  5.  In  all  cases  Eq.  (52)  reproduces  the  original 
values  to  with  10%  and  usually  much  better. 

In  the  past  classical  formulas  have  been  used  to  calculate  the  ionization  cross  sections 
from  which  one  can  obtain  the  three-body  recombination  constants  and  the  ionization  rate 
coefficients  [2].  For  instance  the  classical  formula  of  Thomson  gives  the  ionization  cross 
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IONIZATION  RATE  CONSTANTS 
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THREE-BODY  RECOMBINATION  CONSTANTS 
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Table  5.  Values  of  a  and  t]  in  K{c  —*  p\T)  =  oT  for  two  different 
temperature  regions. 


State 

T  =  5x103  -  2x10“  “K 

T  =  2x10“ 

-  2x10®  “K 

a 

V 

a 

V 

3s 

1.134[-24]“ 

1.044  999 

1.806[-22]“ 

1.553  714 

3p 

5.446[-24] 

1.025  274 

1.015[-21] 

1.549  876 

3d 

3.406[-23] 

1.064  035 

6.800[-20] 

1.596  219 

4s 

3.328[-23] 

1.243  no 

2.538[-21] 

1.679  510 

4p 

1.642[-22] 

1.242  221 

1.516[-20] 

1.698  688 

4d 

6.520[-22] 

1.271  318 

5.305[-20i 

1.715  270 

4f 

6.287[-22] 

1.227  865 

9.058[-20j 

1.729  279 

5s 

2.250[-22] 

1.348  549 

1.193[-20] 

1.748  639 

5p 

9.654[-22] 

1.351  053 

4.800[-20] 

1.745  048 

5d 

2.663[-21] 

1.360  183 

1.385[-19] 

1.759  402 

5f 

3.134[-21] 

1.337  922 

2.304[-19] 

1.771  964 

5g 

3.176[-21] 

1.310  384 

3.499[-19] 

1.785  248 

6s 

7.302[-23] 

1.470  302 

3.356[-20] 

1.793  332 

6p 

2.315[-21] 

1.387  725 

1.236{-19] 

1.789  453 

6d 

5.457[-21] 

1.392  790 

3.057(-19] 

1.779  628 

6f 

7.718[-21] 

1.392  499 

4.534[-19] 

1.804  076 

6g 

9.429[-21] 

1.382  759 

6.855[-19] 

1.815  851 

6h 

1.090[-20] 

1.374  620 

9.474[-19] 

1.825  962 

“The  number  inside  the  bracket  indicates  the  power  of  10. 
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section  (t(E)  at  incident  electron  E  as 

<t{E)  =  {ne*)  -  E-^)  ,  (53) 

where  Ep  is  the  ionization  energy  of  state  “p”.  For  comparison  we  have  computed  the  three- 
body  recombination  constants  by  means  of  this  classical  empirical  procedure.  In  applying 
Eq.  (53)  to  calculate  the  ionization  cross  sections,  we  use  the  experimental  value  of  Ep  rather 
th2in  the  hydrogenic  value.  The  resulting  classical  values  for  the  three-body  rea>mbination 
constants  axe  larger  than  the  quantum-mechanical  counterparts  at  low  electron  temperature. 
For  example  at  5000“K  the  classical  values  are  larger  by  as  much  as  20%.  Conversely  at 
100,000®K,  the  classical  values  are  smaller  than  quantum-mechanical  ones  by  as  much  as 
35%. 
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PART  3 


ELECTRONIC  STRUCTURE  OF  THE  HIGHLY  EXCITED  STATES 
OF  THE  OXYGEN  MOLECULE 

Rydberg-state  wave  functions  that  converge  to  the  OJ  ground  state  have  been  calculated 
for  five  difference  values  of  the  intemuclear  distance  R  (1.0,  1.117,  1.207,  1.3,  and  1 .1  ,d) 
for  the  total  molecular  symmetry  of  OJ (X*n^)(n^)'’^nj,u  for  (n£)  from  (3,0)  to  (10,0)  Ijy 
a  Hartree-Fock  scheme  in  which  the  Coulomb  and  exchange  interaction  of  the  Rydberg 
electron  with  the  OJ  core  is  treated  rigorously.  The  wave  fimction  of  the  OJ (X^Ilg)  ion  core 
is  obtained  by  a  multiconfiguration  self-consistent  field  method.  Comparison  of  the  vertical 
excitation  energy  with  the  available  experimental  data  shows  a  very  good  agreement  for  the 
lowest  npffu  ^IIu  state,  but  a  ten  percent  discrepancy  is  found  for  the  Zsag  states.  The 
ionization  energies  are  found  to  be  insenstive  to  the  core  wave  functions.  Variation  of 
the  ionization  energies  with  the  intemuclear  distance  R  is  generally  a  few  percent  or  less 
from  R  =  1.0  to  1.4A. 

Next  we  study  the  O2  Rydberg  states  that  converge  to  the  excited  states  of  the  parent 
O2  ion,  i.e.,  the  and  c^E~  states  of  O2.  The  wave  fimctions  of  the  Rydberg 

electron  are  again  determined  by  a  Hartree-Fock  scheme.  The  calculations  cover  (nf)  from 
(3.0)  to  (10,7)  with  m=0  and  1.  For  the  intemuclear  distance  we  choose  the  equilibrium 
distances  of  the  parent  excited  O2  states  and  also  R  =  1.2A.  The  results  are  compared  with 
photoabsorption  and  electron-impact  energy-loss  experiments.  Good  agreement  is  seen  for 
the  nsffg  np^u,  and  npTTu  series  but  the  nAag  states  show  larger  discrepancy. 

Details  of  our  works  have  been  published  in 

(a)  Journal  of  Physics  B:  Atomic,  Molecular  and  Optical  Physics, 

Volume  21,  pp.  1155-1165  (1988), 

(b)  Physical  Review  A,  Volume  39,  pp.2367-2373  (1989). 
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PART  4 


ELECTRON-IMPACT  EXCITATION  CROSS  SECTIONS  OF  THE 
RYDBERG  OXYGEN  MOLECULES 

We  have  formulated  a  scheme  for  calculating  electron-impact  excitation  cross  sections 
for  the  Rydberg  states  of  the  oxygen  molecule  by  means  of  the  method  of  close  coupling. 
This  was  based  on  our  earlier  works  on  the  application  of  the  method  of  close  coupling  to 
calculate  cross  sections  for  electron-impact  excitation  of  electronic  states  of  the  H2  and  N2 
molecules.  In  the  earlier  works  we  were  dealing  with  excitation  into  the  low-lying  electronic 
states  of  N2  and  H2,  the  wave  functions  of  which  were  expressed  as  linear  combinations  of 
atomic-like  orbitals  at  e2w:h  nucleus,  i.e.,  in  the  two-center  form.  In  the  present  work  we 
consider  excitation  into  the  Rydberg  states  of  O2,  and  their  wave  functions  are  expressed 
in  the  one-center  form.  A  considerable  modification  had  to  be  made  in  order  to  adapt  the 
final-state  wave  functions  to  the  one-center  form.  We  have  used  this  method  of  close  coupling 
to  calculate  the  cross  sections  for  exciting  O2  from  the  ground  electronic  state  (X®E~)  to 
the  [Oj  (b^Ej^)](3p<ru)^Eu  Rydberg  state  at  60  and  80  eV  and  the  results  agree  with  those 
of  the  Bom  approximation  to  within  4%.  This  indicates  that  the  Bom  approximation  is 
sufficiently  accurate  for  this  work. 

We  have  made  a  comprehensive  study  of  electron-impswt  cross  sections  for  exciting  the 
oxygen  molecule  from  the  ground  electronic  state  to  the  nsag,  npa^,  ndcTg,  nfcTu,  npTTu, 
ndTTg,  and  nf7r„  Rydberg  states  associated  with  the  and  02(b'*E")  cores.  The 

cross  sections  are  calculated  by  the  Bom-Ochkur  and  Bom-Rudge  approximations.  The 
wave  function  of  the  X^E“  groimd  state  of  O2  is  determined  by  a  multiconfiguration  self- 
consistent-field  calculation.  For  each  Rydberg  orbital  associated  with  the  O2  )  core, 

we  have  a  singlet  and  a  triplet  state,  whereas  each  Rydberg  orbital  associated  with  the 
core  gives  rise  to  a  triplet  and  a  quintet  state.  Since  the  ground  sate  of  the  O2 
molecule  is  a  triplet,  we  have  both  spin-conserving  and  spin-changing  excitation.  At  high 
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energies  the  excitation  function  for  the  spin-conserving  excitation  are  found  to  conform  to 
the  Bethe-tjrpe  asymptotic  forms,  whereras  the  excitation  functions  for  the  spin-changing 
excitation  show  the  characteristic  E“^-dependence  at  large  E.  Once  the  cross  sections  are 
available  up  to  200  eV,  one  can  easily  extrapolate  to  higher  energies. 

Within  each  Rydberg  series  of  a  given  I,  m,  and  Oj  core  state,  the  shape  of  the  excitation 
function  is  nearly  identical  for  the  various  n-member  except  the  first  member  (the  lowest 
level)  for  which  the  excitation  function  has  a  somewhat  different  shape  with  the  peak  at  a 
slightly  lower  energy.  Furthermore,  for  each  series,  cross  sections  of  the  various  n-niernbers 
(except  the  first  member)  exhibit  an  n  dependence  of  the  form  n”“  with  a  varying  from  one 
series  to  another  in  the  range  of  1.47  to  4.60  in  contrast  to  the  more  universal  n”^  dependence 
for  the  hydrogen  atom. 

Details  of  this  work  have  been  published  in  the  Physical  Review  A,  Volume  42,  pp. 
4391-4394  (1990). 
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PART  5 


EMISSION  CROSS  SECTIONS  FOR  THE  SECOND  NEGATIVE  BAND 
SYSTEM  OF  OXYGEN  PRODUCED  BY  ELECTRON  IMPACT 

We  have  meafured  the  absolute  optical  emission  cross  sections  for  the  O2  — ► 

X^Wg^v")  second  negative  band  system  for  v’  from  0  to  13  produced  by  electron  impact 
with  O2  in  the  ground  state  =  o)  for  incident  electron  energy  from  0  to  hOO  t-V 

A  monoenergetic,  narrow  electron  beam  of  current  about  90/iA  is  passed  through  a  collision 
chamber  filled  with  oxygen  gas  at  a  presure  about  10  mTorr.  By  measuring  the  absolute 
intensity  of  the  emission  for  the  various  v'  — >  v"  bands,  we  have  determined  the  absolute 
optical  emission  cross  sections  for  some  28  bands.  There  are  a  number  of  emission  bands 
that  are  too  weak  to  be  detected  in  an  electron-beam  experiment.  For  some  of  these  weak 
bands,  we  used  an  indirect  method  to  obtain  the  emission  cross  sections.  In  ein  oxygen-gas 
discharge  the  second  negative  emission  bands  are  much  stronger,  but  the  OJ  ( A^Ilu,  v')  levels 
are  populated  by  mechanisms  that  are  more  complicated  tham  direct  electron  excitation. 
However,  regardless  of  the  mechanisms  of  population,  the  ratio  of  the  intensity  of  two  bands 
originating  from  the  same  upper  state  is  equal  to  the  ratio  of  the  corresponding  optical 
emission  cross  sections,  i.e., 

^  (»■'  -  O  ^  (”'  -  <)  ,54, 

/  (v'  — >  «(,')  (v'  — >  V)l) 

Suppose  we  have  meaisured  the  emission  cross  section  of  a  strong  band,  v'  — >  vH,  but  the 
v'  — >  band  is  too  weaJc  to  detect  in  an  electron-beam  experiment.  Thus  we  resort  to  an 

oxygen  discharge  in  which  all  the  bands  are  much  stronger,  to  determine  the  intensity  ratio 
[l(u'  — ♦  v*J)/I{v*  — >  Combining  this  intensity  ratio  with  the  measured  cross  section 

— >  vj^),  we  can  then  obtain  the  cross  section  for  the  weaker  band  v'  — *  v".  In 
this  manner  we  have  obtained  the  emission  cross  sections  for  a  number  of  weaker  bands. 
From  these  optical  emission  cross  section  data  we  determine  the  apparent  cross  sections  for 
ionization-excitation  to  14  vibrational  levels  of  O3  (i4’n„,  v')  are  determined. 
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We  have  measured  the  variation  of  emission  radition  intensity  of  the  5— ►  4  band  of  the 
Second  Negative  System  with  the  radial  distance  from  the  center  of  the  electron  beam.  The 
result  is  shown  in  Figure  15.  There  is  no  significant  ttiiling  of  the  intensity  with  increasing 
distance.  Therefore,  this  measurement  does  not  show  evidence  of  a  significant  long-life 
component  in  the  population. 

Details  of  our  study  of  the  cross  sections  for  the  second  negative  band  system  have  been 
published  in  Physical  Review  A,  Volume  38,  pp.4537-4545  (1988). 
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Position  (mm) 

Variatioi  of  the  emissiot,  radiation  Ir. tensity  cf  the 
band  of  the  .Second  Negative  System  of  the  oxyget 
molec.jl»  with  the  radial  distance  from  the  center  of  the 
electron  beam. 


PART  6 


FORMATION  OF  HIGHLY  EXCITED  STATES  BY  ELECTRON 
RECOMBINATION  WITH  ATOMIC  AND  MOLECULAR  OXYGEN  IONS 

Cross  sections  for  rauliative  recombination  of  electrons  with  groimd-state  oxygen  ions  to 
form  oxygen  atoms  in  the  0[(^5)n^]^’®L  states  have  been  calculated  for  nl  up  to  (15,1  1) 
at  19  electron  energies  from  0.0345  to  2.585  eV.  The  wave  functions  for  the  bound  states 
and  the  continuum  states  are  determined  by  the  Hartree-Fock  method.  Extrapolation  of  the 
cross  sections  to  n>15  enables  us  to  determine  the  radiative  recombination  rate  coefTicients 
at  different  temperatures.  Analysis  of  the  individual  cross  sections  show  that  our  calculated 
cross  sections  are  very  close  to  those  obtained  by  using  the  hydrogenic  approximation  for 
£  >  2.  This  is  because  the  centrifugal  barrier  keeps  the  active  electron  from  sampling  the 
bound  electrons  of  O'*",  so  that  only  the  net  positive  charge  of  the  O'*'  ion  is  important. 
Therefore,  the  cross  sections  that  we  have  obtained  for  electron  radiative  recombination 
with  O'*"  to  form  oxygen  atoms  in  the  n£  states  with  £  >2  axe  also  applicable  to  radiative 
recombination  of  electrons  with  Oj  to  form  the  corresponding  n£  Ryderg  states  of  the  O2 
molecule.  More  details  of  this  work  can  be  foimd  in  our  paper  which  was  published  in 
Physical  Review  A,  Volume  43,  pp.3433-3439  (1991). 
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